Unified study of glass and jamming rheology in soft particle systems 
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We explore numerically the shear rheology of soft repulsive particles at large volume fraction. 
The interplay between viscous dissipation and thermal motion results in multiple rheological regimes 
encompassing Newtonian, shear-thinning and yield stress regimes near the 'colloidal' glass transition 
when thermal fluctuations are important, crossing over to qualitatively similar regimes near the 
'jamming' transition when dissipation dominates. In the crossover regime, glass and jamming sectors 
coexist and give complex flow curves. Although glass and jamming limits are characterized by similar 
macroscopic flow curves, we show that they occur over distinct time and stress scales and correspond 
to distinct microscopic dynamics. We propose a simple rheological model describing the glass to 
jamming crossover in the flow curves, and discuss the experimental implications of our results. 

PACS numbers: 62.20.-x, 83.60.La. 83.80.Iz 



The emergence of solidity in disordered assemblies 
of repulsive particles is a well-known phenomenon [l| 
which remains difficult to understand at a fundamental 
level [2|, |3[. When compressed, a colloidal suspension 
undergoes a 'glass transition' from (metastable) thermal 
equilibrium [^, as observed experimentally for a broad 
spectrum of particle types 15| . For colloidal hard spheres 
suspended in a solvent of viscosity rjs , the shear viscos- 
ity, rjT, is a universal function of the packing fraction 
(y9; rjx/rjg = G{(p), independently of e.g. particle size [6[- 
Solidity also emerges far from equilibrium upon com- 
pressing non-Brownian suspensions of repulsive particles 
across the 'jamming transition', as in foams or granular 
materials [3|, |5| . The viscosity 770 of a non-Brownian hard 
sphere suspension is again universal, 77o/'7s = J{f) [l|- 
Depending on the community and the particular system 
at hand, rheologists use a broad variety of functional 
forms and empirical models for G{ip) and J{^), while 
underlying physical processes for both limits are often 
not distinguished in rheology textbooks [8|. Our aim is 
to determine if and how these two ideal limits are interre- 
lated, addressing also the non-linear rheological regimes 
and the additional effects of particle size and softness. 

Glass and janrming transitions share important sinr- 
ilarities, in particular at the rheological level. In both 
cases, solidity emerges near a 'critical' volume fraction 
below which the material is a fluid whose viscosity in- 
creases rapidly with (p. The anrorphous solid at large 
density responds clastically for small deformation, but 
flows when a stress larger than a yield stress is applied [l[ . 
The dynamics becomes very heterogeneous near the criti- 
cal density; its spatial correlations are usually interpreted 
by appealing to underlying phase transitions [9|, though 
the nature of these remains a subject of debate [3, l3|. 
Based on these sinrilarities, a unified jamming phase di- 
agram has been proposed where thermal and athermal 
systems appear as a single 'jammed' phase lCllll|. 

At the theoretical level, recent results have clarified the 



relation between glass and jammed phases |12h14|. sug- 
gesting that the jamming transition occurs well inside the 
non-ergodic glassy phase. For systems of soft repulsive 
particles, as studied below, a glass transition line Tg{ip) 
separates the fluid (at high T, low ip) from the glass 
(at low T, large (p), while the jamming transition occurs 
upon compression at T = inside the glass phase [14| . 
However, these static calculations shed little light on ei- 
ther dynamical properties or rheology [l5| . Although this 
theoretical scenario appears in broad agreement with nu- 
merical work 



11^, 



16l |. glass and jamming transitions are 
typically located using different sets of methods and ob- 
servables, normally requiring extrapolation |16j . Similar 
ambiguities exist in experimental work where e.g. esti- 
mates for the location of the colloidal glass transition 
cover the range ipc ~ 0.57. . .0.635, depending on how 
the divergence of G{ip) is extrapolated [J, [g, [l3|- In the 
same vein, data for the divergence of J{ip) for athermal 
suspensions he in the range tpj « 0.585 . . . 0.66 ITI- Ildllq . 

In this paper, we argue that a clearer picture emerges 
when the non-linear rheology of both thermal and ather- 
mal suspensions is considered. We use computer simu- 
lations to investigate the flow properties of concentrated 
assenrblies of soft repulsive particles, and vary the relative 
strength of thermal fluctuations and viscous dissipation 
to study the crossover from thermalized suspensions (rel- 
evant to soft colloids) to purely athermal ones (relevant to 
jammed solids) within a single computational framework. 
This setting allows us to observe and disentangle multi- 
ple rheological regimes within a single system, establish- 
ing in particular unambiguously that the increases of the 
shear viscosities »7t(<p) and rio{(p) upon compression are 
unrelated. This has inrportant consequences for the janr- 
ming phase diagram of soft particles. Although glass and 
jamming limits are characterized by similar macroscopic 
flow curves, we also show that they in fact occur over 
well-separated time and stress scales and correspond to 
qualitatively different microscopic dynamics. 



We analyze theoretically the behaviour of sheared as- 
seniblies of soft repulsive particles immersed in a solvent, 
such as star polymers, microgels, or dense emulsions Q. 
The simplest way to model these systems at large pack- 
ing fraction is to ignore hydrodynamic interactions and 
consider only pairwise repulsion between particles, such 
as V{r) = e(l — r/a)"6(a — r), where &{x) is the Heavi- 
side function, and a is the particle diameter [19| . We use 
molecular dynamics simulations to study the steady state 
rhcology of harmonic spheres, a = 2, in a simple shear 
flow. We simulate the following Langcvin dynamics. 
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where ri and yi represent the position and the y- 
coordinate of particle i, respectively, and e^ is the 
unit vector along the x-axis. The damping coefB- 
cient, ^, which accounts for viscous dissipation, and 
the random force fi describing thermal fluctuations 
obey the fluctuation-dissipation relation, {fi{t)fj{t')'^) = 
2kBT^dijl6{t — i'), where fc^ is Boltzmann's constant. 

The Langevin dynamics in Eq. ([1]) is characterized by 
two microscopic timescales: tq = ^a^/e controls the dis- 
sipation, while TD = ^0? /{ksT) ~ Toe/{kBT) sets the 
timcscale for Brownian motion. Therefore, td and tq 
are comparable at high temperatures but become well- 
separated when fcfiT <§; e, with tjj ^ tq. The shear 
rate 7 introduces a third timescale, 7"^, from which the 
Peclet number is defined, P^. = ^t^. Timescale separa- 
tion at low T allows us to separately explore the ther- 
mal regime at small 7, P^ <C 1, where Brownian mo- 
tion is relevant, and the athermal regime at larger 7, 
Pe S> 1, where the suspension is non-Brownian. In the 
alternative SLLOD dynamics frequently used to shear 
suspensions [20|, inertia is included and the thermostat 
is the only source of dissipation. Thus, the only acces- 
sible regime is Po < 1 (with now tu = \/|p%)> and the 
T — !■ limit is unphysical. By contrast, Eq. ([T|) can be 
simulated at T = (i.e., Po = 00), where the dynamics 
becomes similar to earlier studies of the jamming transi- 
Thus, Eq. (H)) allows us to study thermal and 
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tion 

athermal systems under shear in a unified manner. 

We study a Zd system of TV = 10'^ harmonic spheres, 
using a 50:50 mixture of spheres with diameter ratio 1.4 
to avoid crystallisation [ly, |21|. We measure length in 



units of the small particle diameter, a, time in units of tq 
and temperature in units of e/fcs- We integrate Eq. ([T]) at 
constant 7 using Heun's method with Lees-Edwards pe- 
riodic boundary conditions 120| over a typical simulation 
time of at least IO/7. We evaluate the xy-component of 
the shear stress, a^ using the Irving-Kirkwood formula, 
and deduce the shear viscosity, t] = aj^. The stress and 
viscosity units are respectively e/a^ and ^/a. For thermal 
simulations at low temperatures, we combine data from 
Langcvin and SLLOD dynamics to broaden the range 
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FIG. 1: Flow curves at different temperatures and densities, 
shown as (7(7) (left) or 77(7) (right). Diamonds (left) mark 
the state points analyzed in Fig. [S] Circles (right) indicate 
thermal (closed) and athermal (open) viscosities reported in 
Fig. [2I Flow curves are shown in blue (dot-dashed) when 
thermal Newtonian behaviour is observed, in red (dashed) 
when thermal effects are negligible, in black otherwise. 



of shear rates towards very low Peclet numbers, where 
SLLOD is more efhcient. We have checked that both 



methods yield equivalent results at equal Pg values [31 



We first characterize the macroscopic flow properties 
at T = 10~*. Since ty) = lO"', the data in Fig. [IJi mainly 
cover the thermal range, Po < 1. Accordingly, the re- 
sulting flow curves are typical of dense fluids sheared 
across the glass transition [15|, |22|, [23 1. Briefly, for den- 
sities ip < ipG ~ 0.61, flow curves are Newtonian at low 
shear rates, while for larger 7, the external flow acceler- 
ates structural relaxation leading to shear-thinning. The 
Newtonian viscosity, 77t(v), increases rapidly upon com- 
pression towards ipc- Above ipc, the linear viscosity is 
too large to be measured, and the system behaves instead 
as a solid (a glass) with a finite yield stress, defined as 
(7y = lini^_>.o (7(7). Both the shear viscosity riT{ip) and 
the yield stress ay if) can be used to locate the glass 
transition, which corresponds to either the divergence of 
rjT, or the emergence of a finite cry, see Fig. [2] 

At significantly lower temperature, T — 10^^, the 
above phenomenology persists as long as Po < 1, which 



now corresponds to very low shear rates, 7 < t 
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in Fig.[Tl3. Thus, we can determine a Newtonian viscosity 
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FIG. 2: (a) Shear viscosities r]T and ?7o and their distinct hard 
sphere limits G((/?) [l7| and J{}p) [23|. (b) Density dependence 
of yield stress for different temperatures, including the T = 
limit, (c) Same data as in (b) in a 'glass-jamming phase 
diagram'. 



r]T{sp) for 7 — !> and ^p < ipa ~ 0.59, and a finite yield 
stress above ipc ■ Note that ipc decreases slowly with de- 
creasing T, see Ref. [16[. Because td is very large, there 
now exists a broad 7 window where Po 3> 1 and ther- 
mal fluctuations play little role. Surprisingly, the data in 
Fig. [Ud show that for a range of densities above the glass 
transition, 0.59 < tp < 0.63, the system flows as a simple 
Newtonian fluid when Pe ^ 1- This allows us to define 
a second viscosity, 770, that grows upon compression to- 
wards Lpj ~ 0.64. Finally, for Lp > Lpj, the flow curves 
are mainly characterized by a yield stress, (Ty((p). The 
shear viscosities i]t{^) and r]o{(p) (Fig. [2^) obey clearly 
distinct behaviours. While the growth of rjT reflects the 
approach to the glass transition, ipc ~ 0.59, the increase 
in r]o is separately controlled by the jamming transition, 
(fj ~ 0.64. Given that both viscosities are deflned over 
distinct density and shear rate regimes, and can be simul- 
taneously observed at this temperature, it is clear that 
they reflect distinct phenomena, even without extrapo- 
lation to locate ipG and ifj more precisely. Correspond- 
ingly, the evolution of cryiv) in Fig. [2b is influenced by 
both transitions, since solidity emerges near ipQ^ but ay 
increases significantly near ifj > ipc- This is consistent 
with the idea that jamming mainly affects the very low 
temperature properties of the glass phase [14 1 . 

Finally, the rheology at T = corresponds to Po = cxd, 
and so the glass physics cannot operate. Despite this 



ones obtained at T = 10~^ in Fig. [T^. They are charac- 
terized by a Newtonian viscosity rjo at small 7 and low 
density, Lp < Lpj k, 0.64, while a yield stress emerges upon 
compression, Lp > (pj. These data are fully equivalent 
to previous rheological studies of the athermal jamming 
and indeed near that transition can be 
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transition 

collapsed using the same critical scaling. The qualitative 
similarity between flow curves in Figs. [TJi and[lJ: has cre- 
ated confusion in the literature [15[ , where data obtained 
for systems undergoing the glass transition have been in- 



correctly analyzed in the athermal scenario of [21 



The shear viscosities rjx and 770 a-re measured over 7 
windows that become well-separated at low T and Fig. [2^ 
emphasizes that the difference between the two functions 
increases as T decreases, ruling out a smooth convergence 
of rjT to rjQ for T — > 0. Instead, we find that as T — ;> 0, 
riT(Lp) follows the same density dependence as the equi- 
librium relaxation time of the corresponding hard sphere 
fluid IITII. while 770 is well described by an algebraic diver- 
gence [l8|. Our results establish that the functions G((p) 
and J{lp) controlling Newtonian flow in the hard sphere 
limit arc conceptually and quantitatively distinct. 

The yield stress is another highly sensitive indicator 
of the differences between glass and jamming transitions, 
see Fig. ^p. At finite T, solid behaviour emerges near 
LPG ~ 0.59 . . . 0.61, which agrees with equilibrium dynam- 
ics studies [16( . The yield stress then increases smoothly 
with ip up to Lpj « 0.64 where its density dependence 
changes dramatically. Also, while ay scales with T be- 
low Lfij, it is of order unity (in our units of e/a^) above, 
with only a weak T-dependence scaling approximately as 
~ ((^ — ifij). Consistent with this picture, more detailed 
analysis shows that (Ty{ip,T) follows scaling behaviour 
near (pj very similar to the one derived in Ref. |24| for 
the pressure [32|. Thus glass and jammed states, hav- 
ing distinct physical origins, also display distinct stress 
scales, and remain well-separated even as T — > in the 
'glass-jamming phase diagram' shown in Fig. [2]:. Note 
also that while the glass transition occurs at finite T in 
the unsheared system, the jamming transition exists at 
T = only, so that these two limits never coexist. 

The complex flow curves shown in Fig. [T] can be mod- 
elled by assuming that the stress is an additive combina- 
tion, 17(7) = ac + uj + rjaj, of contributions from glass 
and jamming physics, and from the solvent. The sim- 
plest model for the glass contribution incorporating the 
appropriate time and stress scales is 



C^G 



Sg('^) 



(fcsT/aS) [1 + (7Td)^]-i + iiTDfciv)]- 



(2) 



complete change of regime, the corresponding flow curves 
shown in Fig. [TJ; appear qualitatively very similar to the 



with the dimensionless stress Sg(<^) = const. + {p — 
lPg)"'{'-Pj — Lp)~^Q{Lp — lpg)Q{lpj — Lp), and a dimension- 
less timescale, tg((^), diverging at ipc, e.g. as tg((^) ^ 
(lpq — Lp)~'^ . In this model, the viscosity diverges at LpQ, 
t^t/Vs k tg{ip) (with r]s = C/a), and a flnitc yield stress 
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FIG. 3: Mean-square displacement at (/? = 0.62 and 7 = 
10^^ • • • 10^^ and different temperatures. Caged dynamics 
(dashed) is only observed at finite T and Pc <^ 1. 



emerges above, ayi^ > fa) = [ksT / a^)TiQ{ip) which 
diverges near ipj. For the jamming contribution, we use 
a model consistent with the scaling discussed in |21| . 



aj/(e/a3) = tj{^) + [{^t^)-^' + (7T0 J(^))- 



(3) 



where we take Sj(<y9) ^ [ip — fj)" 6(</' — ^j) and 
J (if) '^ (if J — ip)^'^ . The entire set of data shown in 
Fig. [T] can be reproduced nearly quantitatively usin g ex - 
ponents consistent with existing results: a = 0.7 [25j . 
/3 = 0.3 m, ,5 = 0.8 [13 7 = 2.2 ^, /3' = 0.4, a' = 1.2 
and 7' = 2.0 [2l|. Our empirical model, Eqs. dUIS]), em- 



phasizes that glass and jamming physics take place over 
distinct time and stress sectors, and can be addressed by 
independent theoretical means. Differi ng pred ictions for, 
e.g., whether the onset of yield stress [l5|,l26| is continu- 
ous, as in theories of driven athermal systems [27|,|28[, or 
discontinuous as for driven glasses 25| then make physi- 
cal sense. 

The distinct nature of thermal and athermal regimes 
is apparent also in microscopic dynamic correlation func- 
tions. We plot the mean-squared displacements, i?2(i) = 
Unit) - fiiO)\'^) in Fig. El for fixed p = 0.62 and different 
T and 7. In the glass regime, T = 10~^, i?2(i) displays 
short-time diffusion, caged dynamics at intermediate 
times, and shear driven diffusion at long times [2^ |23| . 
At T = 0, we obtain very different, superdiffusive and 
diffusive, behaviours, as discussed in Ref. [29|. For in- 
termediate temperature, T = 10~^, Fig. [3] shows a clear 
crossover between thermal and athermal regimes: while 
caged dynamics is observed for low Pa, superdiffusive 
motion is obtained at large Pc- Therefore, the thermal- 
athcrmal crossover observed in the macroscopic rheology 
in Fig. [T] originates from a similar crossover at the micro- 
scopic level. While macroscopic flow curves in Fig. [T^-c 
can easily be confused, microscopic observations as in 
Fig. [3] provide a clear qualitative distinction. 

We have used temperature to study the crossover be- 
tween two limits, while experimentalists might equiva- 
lently tune particle softness. We note from Figs. \V^ and 
[2}j that for temperatures above T ^ 10^^, correspond- 
ing to thermal particle compressions (1 — r/a) ^ T^/"^ 
of only 10"^ . . . 10"^, the T = physics has little influ- 



ence on the rheological behaviour. This suggests that the 
jamming transition cannot be studied using soft colloids 
unless T / e is extremely small. A second relevant experi- 
mental parameter is the particle size setting the timescale 
for Brownian motion, with td '^ Is for particles of 1/im. 
This implies that the present thermal-athermal crossover 
should be observable in experiments by tuning the parti- 
cle size in the range 1-10/xm. Yield stress data for emul- 
sions [30| seem consistent with the data shown in Fig.[2]3, 
but further studies are needed to confirm our predictions. 

In conclusion, we have used shear rheology to study 
the relationship between glass and jamming transitions. 
While both correspond to the emergence of solid be- 
haviour as signalled by a finite yield stress, we have 
demonstrated that they occur over stress and time win- 
dows that become well-separated at low temperatures in 
dense repulsive systems. The glass-jamming phase dia- 
gram (Fig. ^) has a scale-separated 'wing' between (pc 
and (fj, so that the glass transition line does not extrapo- 
late to the jamming point for T — >■ 0. The two transitions 
can only be observed separately in these two distinct lim- 
its, i.e. on the glass line and at the jamming point. Any 
other state point in the phase diagram is in principle af- 
fected by a combination of both phenomena, in a way 
described by the simple theoretical model of Eqs. (|2l[3|). 
This conceptual clarification should help rationalize both 
experimental data and the scope of different theoretical 
approaches. 
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